Low velocity limits of cold atom clocks 



J. Munoz0 I. LizuainQ and J. G. Mug£@ 

Departamento de Quimica-Fisica, UPV-EHU, Apartado Postal 644, 48080 Bilbao, Spain 

Fundamental low-energy limits to the accuracy of quantum clock and stopwatch models in which 
the clock hand motion is activated by the presence of a particle in a region of space have been 
studied in the past, but their relevance for actual atomic clocks had not been assessed. In this 
work we address the effect of slow atomic quantum motion on Rabi and Ramsey resonance fringe 
patterns, as a perturbation of the results based on classical atomic motion. We find the dependence 
of the fractional error of the corresponding atomic clocks on the atomic velocity and interaction 
parameters. 
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I. INTRODUCTION 

There is a sizable amount of work devoted to "quan- 
tum clocks" in which a clock "hand" variable runs dur- 
ing the passage of a particle throu gh a reg ion of space 
0, 0, S i, U H, Sll ©, M El CI Ell, see an ex- 
tensive review in [lj]. Much of this work is related to 
the "Larmor clock" for an electron crossing a potential 
barrier and the determination of tunneling times [l5|, 
but even the timing of a simple freely-moving particle 
poses intriguing questions [l6l |. In principle, and from a 
classical-mechanical perspective, reading the hand after 
the particle has crossed the region (the particle may also 
be put at time t = within the region) should provide 
a measurement of the time that the particle has spent 
there. The method faces though a low energy limit to 
its accuracy because the interaction between the particle 
and the clock perturbs the particle and thus the mea- 
surement. This limitation was formulated by Peres 0, Q 
as a bound of the time resolution of the clock r (the time 
required to distinguish two orthogonal hand states) , 

tE > h, (l) 

for a given particle energy E. The same result follows by 
imposing that the transmission probability be close to 1 
[lfj | ; it has been also found without reference to any clock 
from Kijowski's time-of-arrival distribution p7| . r being 
in this case the time-of-arrival uncertainty and E the av- 
erage incident energy; and quite interestingly, it arises as 
well as a decoherence condition for consistent histories 
in the analysis of Halliwell and Yearsley of the quantum 
time of arrival fl8| . If, encouraged by the agreement of 
so many different approaches, we apply this inequality to 
the energy of a Cesium atom moving at speeds of con- 
ventional Cs clock standards (between 100 and 300 m/s), 
the predicted resolution is ~ 10~ 14 s, which is consistent 
with accuracies achieved in the 70's and 80's but it is 
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an order of magnitude too large for more recent conven- 
tional Cs clocks [l9[. If, instead, we assume a velocity 
of 5 cm/s, one of the smallest speeds projected so far 
for a Cs atom clock in space [2(|, the resulting minimal 
"resolution" is a worryingly high 0.36 /is, many orders 
of magnitude larger than the stated target accuracy of 
10~ 16 s for that clock. In fact current cold-atom foun- 
tain clocks provide accuracies much better than the ones 
that follow from Eq. (fTJ). Clearly the naive application 
of Eq. |T]) is not justified, but the question remains as 
to what the negative effect is, if any, of very low atomic 
motion on atomic clocks. 

Quite independently of the works on "quantum 
clocks" , time-frequency metrology has experienced a phe- 
nomenal progress. Precisely, one of the main factors for 
recent and future accuracy improvement is laser cooling 
and the concomitant use of low atomic velocities [l|| to 
make flight times large and the resonance curve narrower. 
Therefore there is indeed a need to evaluate fundamental 
low energy limits for atomic clocks [U [2^]. The char- 
acteristic quantum clock process described above (the 
particle activating a handle in a space region), is simi- 
lar to the physical basis of an atomic clock but there are 
some important differences as well: in a simple atomic 
clock with a one-field (Rabi) configuration, an atom in 
the ground state passes trough a field region in resonance 
with a hyperfme atomic transition, which induces the in- 
ternal excitation of the atom. Since the ground or ex- 
cited state populations oscillate periodically with Rabi's 
frequency f2 when an atom at rest is immersed in the 
field, it is natural, except possibly for too low energies, 
to relate the passage time to the final internal oscilla- 
tion phase, which can be read out from the excited pop- 
ulation. In the atomic clock, however, the quantity of 
interest is not really the passage time, even though it 
plays an important role, but the frequency of the tran- 
sition. This is determined from the central peak of the 
resonance curve, whose width is inversely proportional to 
the passage time as we shall see in more detail: an anal- 
ysis where the transition of the atom across the field is 
treated classically, which is valid when the atomic kinetic 
energy E = k 2 h/(2m) is larger than other relevant en- 
ergy scales, E >> HA, Ml/2, (A is the detuning, m the 



2 



mass, and kh/m the velocity), shows that the probability 
of excitation to level 2, when the incident atom is in the 
ground state 1 and the field width is I, is given by 



Pti 1 (A) 



SI 2 



A 2 T 



tt 2 + A 2 



(2) 



where T = lm/(hk) is the (classical) in- field flight-time. 
Three relevant properties of P 12 are the mentioned nar- 
rowing of the main peak when T increases, the peak at 
A = irrespective of the atomic velocity (or equivalcntly 
of T), and the symmetry around A = 0. Thus, velocity 
averaging does not affect the central fringe much, so that 
measurement of Pi 2 can be used to steer the frequency of 
an external oscillator close to the reference, ideally unper- 
turbed atomic frequency luq- This steering will be more 
accurate the narrower the peak, i.e., for slower atoms. 
The clockwork then produces the second by counting a 
predetermined number of oscillations of the external os- 
cillator, ps 9.2 x 10 9 for Cs clocks, so that a shift in the 
peak translates into an error in the determination of the 
second of the order of the fractional frequency error. 1 

Most atomic clocks actually use the Ramsey configu- 
ration with two separated fields rather than one, to im- 
prove performance, but the above three properties of Pi 2 
remain true also for the Ramsey scheme. In the Ramsey 
scheme the free-flight time between fields plays the role 
of T as the parameter to be maximized, within techni- 
cal constraints, for narrowing the resonance. We shall 
first restrict the study to the Rabi configuration for sim- 
plicity, and provide later the results for the Ramsey con- 
figuration, since the manipulations are very similar but 
lengthier. Our main aim is to determine the fractional 
frequency error of the clock at "low" incident velocities: 
this is the (absolute value of the) shift of the peak of 
Pi 2 (A) with respect to the ideal value A = divided by 
the atomic transition frequency. The term "low velocity" 
here is subtle and must be clarified: we shall consider ve- 
locities where the semiclassical result (with the peak at 
A = 0) has to be corrected perturbatively because of 
quantum motion, but such that the overall interference 
pattern keeps essentially the same form as for the semi- 
classical limit. Physically this means that reflection from 
the field regions is still very unlikely and E » fth. This 
will enable us to find the corrections by solving first the 
quantum dynamics exactly and then, taking the domi- 
nant and first correction terms in "large fc" expansions. 
An example of the perturbative regime is provided in Fig. 
[TJi, where the quantum curve resembles, but it is slightly 
perturbed with respect to, the semiclassical curve. We 
insist that this regime may in fact involve very small 
velocities according to other criteria, such as the veloci- 
ties in Fig. [T£i, unrealistically low in any actual atomic 
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1 unless of course it is well understood and predictable so that it 
can be corrected by hand. 



FIG. 1: (Color online) Rabi (one-field) scheme, (a) Perturba- 
tive regime: the maximum probability of the exact and semi- 
classical solutions P-^ and Pfl 1 practically coincide at A = 
for a 7r-pulse. Note however the peak asymmetry and the re- 
lated shift at half height. Data: I — 3 mm and v — 5 fim/s. 
(b) Reflection-dominated regime: the central fringe is nar- 
rower and severely deformed (v = 0.5 /im/s, note the different 
scale) , and starts from a critical detuning. In (a) and (b) The 
Rabi frequency is fixed by the 7r-pulse condition £1 — ttv/1. 
The mass of Cs is used in all figures: m(Cs) = 2.2 x 10~ 25 
Kg. 



clock, but useful for illustrating clearly the main features 
of the quantum curve, in particular the asymmetry, and 
the maximum very near A = 0, at least for the naked 
eye (the maximum is in fact shifted but only slightly as 
discussed in more detail later on). 

Reflection from the field becomes relevant for E < 
£l/H, but then the interference pattern is drastically de- 
formed. For a Ramsey configuration this implies a mul- 
tiple scattering scenario with entirely new Fabry-Perot 
resonances highly sensitive to the velocity [2l| and the 
existence of a threshold of Pi 2 with respect to detuning. 
This very extreme regime, see an example in Fig. [T]d, is 
out the scope of current time-frequency metrology so we 
shall instead concentrate on the perturbative regime. 

The objective is to evaluate Pi 2 (A) quantally for the 
transmitted atoms. We shall use as a simplified model 
the effective ID Hamiltonian 

H = £ - RA|2)(2| + ~n(x)(|l><2| + |2)(1|), (3) 

which is obtained, see [23|, H3] , by assuming a semiclas- 
sical effective traveling wave field, and applying rotat- 
ing wave and dipolar approximations. The effective mi- 
crowave field may be realized by two copropagating lasers 
perpendicular to the atomic motion which induce a Ra- 
man transition [2a, k26|- In this case the effective de- 
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tuning A contains, apart from the difference between ef- 
fective laser frequency (difference between the copropa- 
gating laser frequencies) and transition frequency, the re- 
coil term —k\hj '(2m), k^ being the effective wavenumber 
(difference between the two laser wavenumbers) . Finally, 
£l(x) is the effective Rabi frequency that we consider to 
be constant in the illuminated region(s) and zero outside. 

To find the exact solution of Eq. we start by solving 
the stationary Schrodinger equation (SSE) for each zone 
labeled with an index a: on the left of the field (a = I), in 
the field (a — II), and on the right of the field (a = III): 

H a 4> a — E<p a ; (a = I, II, III), (4) 



Hi) Finally, for a = III we have, as for a = I, 
1 



/2tt 



\ia ni e ikx + b UI e- ik *)\l) 



+ (c„^ + d„ je -^)|2>]._ (12) 
In all cases, a — I, II, III, it is useful to introduce the 
vector of coefficients 



v a (Qq , b a , c a , d a ) , 
where T means "transpose" . 



(13) 



where 4> a is a two-component wavevector 

^ = A*) • 



i) For a = I, 



P 

2m 



fiA|2)(2| U ; 



(5) 



(6) 



reduces to two uncoupled equations for <jy-p and (f>^ . The 
general solution (delta-normalized in k space) is thus 



II. TRANSFER MATRICES AND EXCITATION 
PROBABILITY 



Let us consider a plane wave —^=e lkx , incident from 



1 r ikx 

the left (x — > —oo), with the atoms in the ground state 
|1) so that a j = 1, c 2 =0, b IU = d IU = 0. After the 
interaction "barrier" , the transmitted atom traveling to 
x — > oo may stay in |1) or be excited in |2), so that 



V; = (l,rn, 0,ri 2 ) r , 
v m = (in,0,ii2,0) T , 



(14) 



/2tt 



[(a I e ikx + b I e~ ikx )\l) 



(7) 



+ ( Cl e iqx + d I e~ lqx )\2) 



where q = y/k 2 + 2mA/h and k — ^2mE/h 2 . 
ii) For the interaction region, a — II, 



|l-SA|2)<2| + ^fi(|l)<2| + |2)<l|) 
2m 2 



<l>u = E <f>u- 



By diagonalizing the interaction part of the Hamiltonian 
we obtain its wave vectors and eigenstates |A±), 



|A±> = 



1 

Q 



(9) 



where X± 



and Q' = VA 2 + Q 2 . 
From Eq. (|lip and (J9j) we may write the solution of 
the SSE in a — II as 



1 



[(a II e* k + x +b II e-* k + x )\\ + ) 



+ (c„ e * fc -* + d„e-' ifc -*)|A_)], 
with the mode wavenumbers obeying 

h 2 k 2 



h 2 k\ 

2m 



2m 



±- 



(10) 



(11) 



where , tij are reflection and transmission amplitudes 
for left incidence in internal state i — 1,2 and outgoing 
internal state j = 1,2. 

Let us now introduce for the outer regions, a = I, III, 
a vector containing the amplitudes and their derivatives, 



^\x) 



M (x)v a , 



(15) 



and similarly, inside the field (a = II) 



4f( X ) 



M b (x,<j) 1 )v u , 



(16) 



where the dot represents derivative with respect to x. 
The explicit form of the matrices Mq(x) and M&(x) is 



M (a) 









a iqx 



\ 

e -iqx 





\ 



iqe lqx —iqe lqx 



J 

(17) 
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M b (x) = 



1 

Itt 







e 4fc - 


X 


e - lk - 


X 


2A+e ifc + x 


2A + e _ifc + x 


2A_e* fc - 




2\_e- ik - 




a. 


n 


n 




n . 








ik-e tk - 


X 


—ik-e~ lk - 


X 


2ifc + A + e ife + x 


-2jfc + A + e- ifc + x 


2ik-\_e ik - 




-2ifc_A_e~ ife - 




Q 


n 


n 




n 





(18) 



r 



The matching conditions at x\ and x 2 for the wave func- 
tions and their derivatives can now be written as 



M (xi)v I = Mb(xi)v u , 
M b (x 2 )v II = M (x 2 )v III . 



(19) 
(20) 



Eliminating w 7/ from the system above, we end up with 
a transfer matrix T(xx, x 2 , <$>i) which connects the ampli- 
tude vectors of both sides, 

Vj =T(xi,x 2 )v III , (21) 

and is given by the matrix product 

T(xi,ar 2 ) = M {x 1 )- 1 M b (x 1 )M b (x 2 )- 1 M (x 2 ). (22) 



Using Eqs. (|14I17I18[) we may calculate the transmission 
amplitude for passing from |1) to |2) from the left as 



ti 



T 3 i 



T31T13 — T33T11 



(23) 



The exact quantum excitation probability is thus, gener- 
alizing the semiclassical result @, 



if 2 (A) = ||f 12 | 2 . 



III. QUANTUM MOTION SHIFTS 



(24) 



Two different shifts will be considered, A 



(Af) 

Q 



and 



Ag^: the first one is the shift of the maximum of the 
resonance peak with respect to A = 0; the second one 
is the shift of the average between the resonance curve 
points at half height. As hinted by Fig. [Iji they have 
quite a different behavior, the first one being negligible 
for the naked eye whereas the second one is clearly sig- 
nificant. We shall provide analytical expressions for the 
former, and numerical results for the later. 

We shall first calculate the maximum value of P^(A). 
From 



dP- 



dA 







and Eq. 



the maximum must satisfy 

2 



, , * , q h ( dtu * dt 12 

112112 H — 7T~ t ^ +'"12 — T\ 

m \ dA dA 



= 0. 



(25) 



(26) 



The explicit expression in terms of the basic variables 
v, L,Q is long and cumbersome, since each element Ty 
of the matrix T(xi,:z; 2 ) in Eq. (|23|) results from mul- 
tiplying four matrices 4x4. To simplify we shall first 
eliminate by imposing £1 = 4^ , which is the semiclas- 
sical condition for a 7r-pulse at A = 0, and then expand 
ti2 in powers of A up to second order, 



t 



12, A 



7o(fc)+7i(fc)A-t-72(&)A 2 



where 



7o 
71 

72 



[*12]a=0 

dt 12 - 
dA 



dH 



A=0 
12 



dA 2 



A=0 



The probability P^(A) is approximately 
P 1 Q 2 (A) = ||( 7 o+7iA + 72A 2 )| 
Substituting the expansions in Eq. (|26p . 



(27) 

(28) 
(29) 

(30) 
(31) 



I (70 +7iA + 7 2 A 2 )| + 

2< ? 2 -Re[( 7 i + 2 72 A)( 7 o + 71A + 72A 2 )*] = 0. (32) 
m 

Gathering terms, truncating in first order of A and using 



0o(k) = lalo+q — (7i7o+7o7i*) 

6ii k ) = (7i 7o + 7o7i) + V — (71 7* + 727o + 7o72*) 

m 

the peak must satisfy 

#o(fc) + #i(fc)A + 0(A 2 ) = 0, 

which can be solved retaining the dominant order in k to 
find the explicit shift expression 



A (M) (k) - ^ 
[ ' ~ l6m 2 vP 



sin (2kl) . 



(33) 



where v = kh/m is the velocity, see Fig. [5^, and the 
superscript M stands for "maximum" . 

We may proceed similarly for the Ramsey configura- 
tion and write down a new transfer matrix which is of 
course more complicated but in principle explicit with 
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FIG. 2: Comparison between numerical (dashed lines) and 
analytical values (solid lines), of Aq over the transition fre- 
quency ojq for a Cs atom at different velocities, (a) Single Rabi 
7r-pulse with I = 3 mm, Rabi frequency fixed by 7r-pulse con- 
dition, Q. = nv/l. For the analytical results only upper/lower 
bounds of Eq. (I33p have been plotted. The approximations 
are better for higher velocities, (b) and (c) Ramsey n/2 pulses 
with I = 1.5 mm and separated by L, (TV = L/l), Rabi fre- 
quency fixed by 7r/2-pulse condition, Q = nv/2l. 



the aid of a program of algebraic manipulation. This 
new transfer matrix depends on the previous parameters 
and on L, the length between fields. For "large" k, a 
large ratio TV — L/l, and imposing ir/2 pulses at each 
field, now ft — we find the approximate envelope 



(Ag°) 



2^2 



± 



fl 2 TT 



16m 2 vl 3 N 



(34) 



Note that the Ramsey shifts are considerably smaller 
than Rabi shifts because of the TV in the denominator. 
Compared to the Rabi case, the shift for the Ramsey con- 
figuration oscillates with respect to k with a short period 
7r/L, superimposed to the long period n/l, see Fig. [2]d,c. 
These shifts are in any case rather small and, in addition, 
they oscillate rapidly with respect to v so that they will 
average out because of the velocity spread of the atomic 
cloud. 

There is however a second type of shift which is quite 
relevant in practice: it is defined as the average of the 
resonance curve detunings at half height of the central 
peak, Aq 1 ' 1 ' . We numerically find out that the resonance 
peak is asymmetric and that, at variance with the former 
definition, this shift does not oscillate around zero but it 
tends to a constant at "large" k for fixed I and L. 

As shown e.g. in Fig. QJl, this shift is more significant 
than the former, both because of its larger magnitude and 
because of the absence of oscillation with respect to veloc- 
ity, see Fig. [3] Note that the quantum excitation curve 
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FIG. 3: (Color online) Constant values of the fractional fre- 
quency offset for the Rabi (TV = 0) and Ramsey configurations 
versus incident velocity for different values of I and TV = L/l. 
Exact results. 



tends to the semiclassical curve when compared "verti- 
cally" for increasing v and for a given A, — P-f|' — > 0. 
This is compatible with the fact that, when compared 
"horizontally" at half height, the two curves are sepa- 
rated by a shift which stays constant when increasing v. 
If we represent only the central peak, this constant shift 
is less and less visible due to the peak broadening and 
changing scale when v increases. A back-of-the-envelop 
argument may help to gain some intuitive understanding 
of this phenomenon giving relevant dependences. First 
notice that the wavenumber q is not an even function 
with respect to A and provides asymmetry: the excited 
atom finds a lower potential ground with positive detun- 
ing than with negative detuning and the quantum exci- 
tation probability at half height may be approximated as 
(q/k)P scl , see Eq. (gU), or simply q/{2k), larger than 1/2 
for positive detuning and smaller otherwise. Expanding 
for large fc, and approximating A as half the distance 
to the first zero, the vertical difference between quan- 
tum and semiclassical curves at half height is, for the 
Rabi scheme, w 7r/(4fei). On the other hand the slope, 
again estimated from the zero and for the Rabi case is 
w lm/(Kkir). As this slope is also the ratio between ver- 
tical and horizontal shifts, there results a fc-independent 
horizontal shift proportional to l~ 2 for Rabi excitation 2 , 
or to L~ 2 for Ramsey excitation and large TV. This de- 
pendences are confirmed in Figs. 02 [4] and[5l 



2 The full predicted shift is for 2 / (4ml 2 ) which provides remarkably 
good results given the roughness of the argument. 
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FIG. 4: (Color online) Fractional frequency offset versus / and 
N = L/l for v — 5 cm/s. Exact results 
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FIG. 5: (Color online) Fractional frequency offset at half 
height versus / and N = L/l for v = 5 cm/s. Exact re- 
sults(squares) and fits to an inverse square form (dashed). 



Eq. (1) relies on the reflection of particles from the in- 
teracting field. For an atomic beam crossing a field near 
resonance this means E » fl.fi, which is reminiscent of 
Peres' bound. However l/fi does not play at all the role 
of the time resolution of an atomic clock, which is nowa- 
days even much better than the oscillation period of the 
atomic transition, a modest 10 -10 Hz. Motivated by the 
facts that a naive application of Eq. (1) does not provide 
a valid assessment of actual low energy limits of atomic 
clocks, and that nevertheless, the trend to use lower and 
lower atomic velocities requires this type of information, 
we have investigated resonance shifts in atomic clock con- 
figurations of Rabi and Ramsey type due to quantum 
atomic motion. Explicit expressions are provided for the 
shift of the maximum, Eqs. (|33|) and (|34|) . but this ef- 
fect is quite small compared to the peak asymmetry, and 
moreover it would cancell out because of averaging over 
incident velocities. Instead, when the resonance is de- 
fined by the middle point of the resonance peak at half 
height, the peak asymmetry implies a significant shift, 
several orders of magnitude larger than the shift of the 
maximum, and robust with respect to momentum aver- 
aging. Both shifts are smaller in the Ramsey scheme than 
in the Rabi scheme, but even so the half-height asymme- 
try shift would have to be considered in the error budget 
of accurate Raman-Ramsey clocks. It can be reduced 
by increasing the width of the interaction region between 
field and atom, and, in the Ramsey scheme, by increasing 
the free-flight length. The present results motivate fur- 
ther work to investigate similar effects in standing wave 
rather than traveling wave configurations and for smooth 
field intensities. 

Finally, we would like to comment on atomic clock de- 
signs based on trapped ions and illumination pulses in 
time rather than fixed in space. As the ion is trapped, 
it might seem that such a scheme is free from quantum- 
motion shifts but this is not the case. We have provided 
expressions for the corresponding shifts in [271 ] . 



Acknowledgments 



IV. CONCLUSIONS 

The bound for the time resolution of a quantum clock, 
Eq. (1), suggests limitations to the accuracy of atomic 
clocks at low atomic velocities. One of the derivations of 



We are very grateful to R. Wynands for useful com- 
ments on the manuscript. This work has been supported 
by Ministerio de Education y Ciencia (FIS2006-10268- 
C03-01) and the Basque Country University (UPV-EHU, 
GIU07/40). 



[1] A. Baz', Sov. J. Nucl. Phys. 4, 182 (1967). 

[2] V. Rybachenko, Sov. J. Nucl. Phys. 5, 635 (1967). 

[3] A. Peres, Am. J. Phys. 48, 552 (1980). 

[4] M. Biittiker, Phys. Rev. A 27, 6178 (1983). 

[5] A. Peres, Quantum Theory: Concepts and Methods, 

Kluwer Academic, New York, (1995), Chap. 12. 

[6] C. R. Leavens, Solid State Commun. 86, 781 (1993). 



[7] C. R. Leavens and W. R. McKinnon, Phys. Lett. 194, 12 
(1994). 

[8] C. Bracher, J. Phys. B: At. Mol. Opt. Phys. 30, 2717 
(1997). 

[9] C. R. Leavens and R. Sala Mayato, Ann. Phys. Leipzig 
7, 662 (1998). 

[10] Y. Aharonov, J. Oppenheim, S. Popescu, B. Reznik, and 



7 



W.G. Unruh, Phys. Rev. A 57, 4130 (1998). 
[11] A. Casher, B. Reznik, Phys. Rev. A 62, 042104 (2000). 
[12] D. Alonso, R. Sala Mayato, and J. G. Muga, Phys. Rev. 

A 67, 032105 (2003). 
[13] D. Sokolovski, Phys. Rev. A 76, 042125 (2007). 
[14] R. Sala Mayato, D. Alonso, and I.L. Egusquiza, Lect. 

Notes Phys. 734, 235 (2008). 
[15] M. Biittiker, Lect. Notes Phys. 734, 279 (2008). 
[16] J. Muiioz, D. Seidel, and J. G. Muga, Phys. Rev. A 79, 

012108 (2009). 

[17] A. D. Baute, R. Sala Mayato, J. P. Palao, J. G. Muga, 

and I. Egusquiza, Phys. Rev. A 61, 022118 (2000). 
[18] J. J. Halliwell and J. M. Yearslev, larXiv:0903.1957l 
[19] J. Vanier and C. Audoin, Metrologia 42, S31 (2005). 



[20] Ph. Laurent et al., Eur. Phys. J. D 3, 201 (1998); Ch. Sa- 
lomon et al., C.R. Acad. Sci. Paris Srie IV 2, 1313(2001). 

[21] D. Seidel and J. G. Muga, Eur. Phys. J. D 41, 71 (2007). 

[22] S. V. Mousavi, A. del Campo, I. Lizuain, and J. G. Muga, 
Phys. Rev. A 76, 033607 (2007). 

[23] I. Lizuain, S. V. Mousavi, D. Seidel, J. G. Muga, Phys. 
Rev. A 78, 013633 (2008). 

[24] A. Ruschhaupt and J. G. Muga, Phys. Rev. A 76, 013619, 
(2007). 

[25] J. E. Thomas et al., Phys. Rev. Lett. 48, 867 (1982). 
[26] J. Vanier, Appl. Phys. B 81, 421 (2005). 
[27] I. Lizuain, J. G. Muga and J. Eschner, Phys. Rev. A 76, 
033808 (2007). 



